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We investigate a two-dimensional problem of an isolated self-interacting end- 
grafted polymer, pulled by one end. In the thermodynamic limit, we find that 
the model has only two different phases, namely a collapsed phase and a stretched 
phase. We show that the phase diagram obtained by Kumar at al. [Phys. Rev. 
Lett. 98, 128101 (2007)] for small systems, where differences between various statis- 
tical ensembles play an important role, differ from the phase diagram obtained here 
in the thermodynamic limit. 



I. INTRODUCTION 

The physics of single polymer chains in a poor solvent is still not very well understood. 
Away from the 6'-temperature, we know that a polymer will be in either a collapsed or a 
swollen state The mean-square radius of gyration {R^)g scales with chain length 
as (-R^)g ~ const, x A^^*^, where i/ is a critical exponent. At low temperatures, when the 
polymer is in the collapsed state, u = 1/d while at high temperatures an "extended" or 
"swollen coil" state exists where u = 1, 3/4, 0.588 1/2 for d = 1,2,3,4 |l| respectively. 
These values are believed to be exact for d = 1,2 and, with a logarithmic correction, for 
c? = 4. At high temperatures stretching a polymer should produce a state where u = 1 which 
we shall refer to as the "stretched- state. Although there are many theoretical BBQ as 
well as experimental p works on pulling of a collapsed chain, it seems that some issues 
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remain to be fully understood. 

Recently Kumar et al. presented results from exact enumeration of force-induced 
polymer unfolding in two dimension in the context of modeling single molecule experiments. 
For finite systems, they proposed a phase diagram, which has three phases namely a collapsed 
phase, an extended phase and a stretched phase (in addition there is a swollen phase which 
only occurs at F = above the ^-temperature). They found a transition line between 
the stretched state and the extended phase. The phase diagram proposed is presented in 
Fig. [H The lower phase boundary was obtained in both the constant force and constant 
temperature ensembles, and indicates a phase-transition line where polymer goes from the 
collapsed phase to the extended state. However, the upper phase boundary was been only in 
the constant force ensemble and it was proposed that this represents a transition line where 
the polymer goes from the stretched state to the extended state. 

At this point it is pertinent to mention that very few attempts have been made to perform 
single molecule experiments in the constant force ensemble. Only recently did Danilowicz 
et al. [3] perform stretching experiments on single stranded DNA in the constant force 
ensemble. It was observed that at low force, the extension increases with temperature, 

n 

while at high force the extension decreases with temperature. Kumar and Mishra [8| found 
that this decrease is an entropic effect and showed that the upper line is not a true phase 
transition but a crossover effect. 
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FIG. 1: Schematic phase diagram proposed by Kumar et al. [6|]. 
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In this paper we focus our attention on the true nature of the phase-diagram for the 
model in the thermodynamic hmit. We present some further studies of the series data 
trying to gauge the scahng behavior of the model at different points in the phase-diagram. 
While somewhat inconclusive our analysis does indicate that the true phase-diagram (for 
non-zero force) has only two distinct phases for non-zero forces and not three as originally 
conjectured. The extended phase does not exist for non-zero forces and the upper phase 
boundary is a finite-size effect only present when the model is studied at fixed force with a 
variable temperature. 

Hence to really delineate the phase diagram we have also performed Monte Carlo simu- 
lations using the FlatPERM algorithm 91]. We investigate several hypothetical phase dia- 
grams. In particular, we consider the possible scenario that the phases seen are two types of 
stretched phase; one where the polymer is maximally stretched in a rod-like conformation 
and the other where u = 1 though the polymer is not maximally stretched. Using the sim- 
ulation results we are able to confidently deduce that there is no evidence of any additional 
phase or phase transition. 

n 

We would like to emphasise that the "phase diagram" obtained by Kumar at al. [6] for 
small systems may still be relevant in the context of experiments on bio-polymers. In real 
systems of finite size differences between various statistical ensembles do play an important 
role as evidenced not only by this previous study but also by recent experimental work [3]. 
We thus see our discovery of a discrepancy between the finite size "phase diagram" and the 
true infinite size phase diagram as an important contribution to a better understanding of the 
types of finite-size effects that may be of importance to the interpretation and understanding 
of experimental results on small systems. 

In section [TT] we define the model. In section IIIII we first briefly review the evidence 
presented using series analysis to support the conjectured phase-diagram 0, Q and then 
present further results from a more thorough and extensive analysis of the series data casting 
doubt on the upper phase boundary of the proposed phase-diagram. In section[IV]we present 
the conclusive results of the Monte Carlo simulations which do not support the existence of 
any additional phase transitions: we carefully consider various possible scenarios. Finally, 
in section |V] we summarize our final conclusions. 
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II. MODEL 



We model the polymer chains as interacting self- avoiding walks (ISAWs) on the square 
lattice as shown in Fig. [21 Interactions are introduced between non-bonded nearest neighbor 
monomers. In our model one end of the polymer is attached to an impenetrable neutral 
surface (there are no interactions with this surface) while the polymer is being pulled from 
the other end with a force acting in the direction perpendicular to the surface. Note that the 
IS AW does not extend beyond either end-point so the y-coordinate Uj of the j'th monomer 
is restricted hj = yo < yj < = h. 



FIG. 2: The model of a polymer on the two dimensional square lattice pulled by the last monomer. 
The arrow indicates the direction of the pulling force. The dark (red) filled circles on lattice sites 
denote monomers interacting via nearest-neighbor interactions. 

We introduce Boltzmann weights u = exp(— £:/A;bT) and u = exp{F/kBT) conjugate to 
the nearest neighbor interactions and force, respectively, where e is the interaction energy, 
/cb is Boltzmann's constant, T the temperature and F the applied force. In the rest of this 
study we set e = —1 and = 1. We study the finite-length partition functions 



where CN,m,h is the number of ISAWs of length N having m nearest neighbor contacts and 
whose end-point is a distance h from the surface. 
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FIG. 3: The fluctuations in the number of contacts as a function of temperature for fixed force 
F = 0.0 (upper left panel) and F = 0.5 (lower left panel). Each panel contains curves for ISAWs of 
length (from bottom to top) N = 20, 25, . . . , 55. In the upper right panel we show a log-log plot of 
the growth in the peak value of the fluctuation curve with chain length N. The lower right panel 
shows the peak position (critical temperature value) vs 

III. SERIES ANALYSIS 
A. Fluctuation curves and the conjectured phase diagram 



To begin let us recall the type of analysis presented by Kumar et al. p, Il0|. At low- 
temperature and force the polymer chain is in the collapsed state and as the temperature 
is increased (at fixed force) the polymer chain undergoes a phase transition to an extended 
state. The value of the transition temperature (for a fixed value of the force) can be ob- 
tained from the fiuctuations in the number of non-bonded nearest neighbor contacts. The 
fiuctuations are defined as x = {to?) — (^)^, with the fc'th moment given by 
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FIG. 4: The finite-size 'phase diagram' for flexible chains as obtained from the position of the 
peak in the contact fluctuation curves for N = 55. The solid black curve and the dashed curve are 
obtained by fixing the force and varying the temperature. 

In the panels of Fig. [3] we show the emergence of peaks in the fluctuation curves with 
increasing at fixed force F = 0.0 and F = 0.5. In the top right panel we show the 
growth in the peak value as is increased. Since this is a log-log plot we see that the 
peak values grows as a power-law with increasing A^; this divergence is the hall-mark of 
a phase transition. In the lower right panel we have plotted the position of the peak (or 
transition temperature) as a functions of Clearly the transition temperature appears 

to converge to a finite (non-zero) value but the data exhibits clear curvature which makes 
an extrapolation to infinite length difficult. 

In Fig. m we show the force-temperature 'phase diagram' for flexible chains as obtained 
from the peak positions for the finite chains. However, true phase diagram should be ob- 
tained by extrapolating the data to the — oo limit). In Fig. H] we have shown the 
transitions as obtained by fixing the force (black curves). One of the most notable feature 
of the phase-diagram is the re-entrant behavior but this has been studied and explained in 



previous papers 



a, [lOj. The other notable feature is that in the fixed force case we see an 
apparent new transition line from the extended state to the fully stretched state which is 
solely induced by the applied force (the dashed line in Fig. Hj). 
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FIG. 5: Same as in Fig. [3] but for fixed force F = 1.1 (upper left panel) and F = 1.2 (lower left 
panel). 

B. Further series analysis results 

In Fig. [5] we have plotted the fluctuation curves for force F = 1.1 and F = 1.2. The 
curves for F = 1.1 (including the plot of the peak height) looks very similar to the plots 
(see Fig. [6]) for low values of the force. For force F = 1.2 the peak is not very pronounced 
and we are hesitant to even call it a peak. Also when we look at the peak height vs. it 
appears that the curve has two different behaviors for small and large N, respectively. This 
could be a sign of a cross-over behavior. 

One can also study the same transition phenomenon by fixing the temperature and vary- 
ing the force. In the panels of Fig. [6] we show the emergence of peaks in the fluctuation 
curves with increasing N at fixed temperature T = 1.0 and T = 0.5. Again we observe the 
power-law divergence of the peak-value. The only other note-worthy feature is that in the 
plots of the peak position (critical force value) we observe not only strong curvature but we 
actually see a turning point in the curves as A^ is increased. This feature would make it im- 
possible (given the currently available chain lengths) to extrapolate this data. However, we 
do not observe the upper transition line in this study where we have fixed the temperature 
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FIG. 6: The fluctuations in the number of contacts as a function of force for fixed temperature 
T = 1.0 (upper left panel) and T = 0.5 (lower left panel). Each panel contains curves for ISAWs of 
length (from bottom to top) = 20, 25, . . . , 55. In the upper right panel we show a log-log plot of 
the growth in the peak value of the fiuctuation curve with chain length A^. The lower right panel 
shows the peak position (critical force value) vs 1/A^. 

and varied the force. Indeed this is clear from Fig. [6] where at fixed T = 0.5 and 1.0 we see 
only a single peak (giving us points on the red curve in the 'phase diagram' Fig. H]). 

In Fig. [6] the value of the force extends up to F = 2.0 and the upper transition (dashed 
line in the phase diagram) should appear (if present) as a second peak in the fluctuation 
curves of Fig. O The absence of any evidence of a second peak is what leads us conclude 
that we do not see this second transition in the fixed T varying F study. 

In Fig. [7] we have plotted the average extension per monomer {h)/N as a function of 
temperature. In the case of a fixed force F = 1.2 (lower right panel) we note that curves 
for different values of N more or less coincide showing that the average extension scales 
like A^ for all temperatures. We contend that this observed behavior shows that the upper 
joundary is a crossover effect supporting the finding reported recently by Kumar and Mishra 
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FIG. 8: The average number of contacts per monomer {m)/N as a function of temperature T for 
different values of the force. Each panel contain four curves for, from bottom to top, A'' = 25, 35, 
45, and 55. 
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In Fig. [8] we have plotted the average number of contacts per monomer {m)/N as a 
function of temperature. In the case of a fixed force F = 1.2 (lower right panel) we note 
that curves for different values of N more of less coincide showing that the average extension 
scales like N for all temperatures. 



IV. SIMULATION RESULTS 



Our more detailed analysis of the series indicates that the upper phase boundary is 
not a real phase transition.. To further investigate the model we have turned to Monte 
Carlo simulations that allow analysis of longer polymer chains. We have chosen to use 
the FlatPERM algorithm to simulate the model. One advantage of FlatPERM, as a 
"fiat histogram" technique, is the ability to sample the density of states uniformly with 
respect to a chosen parametrisation, so that the whole parameter range is accessible from 
one simulation. This allows us to "see" the phase diagram from one set of results. The cost 
of this however is that the chain lengths that can be simulated accurately are still fairly 
modest. We have performed "whole phase space" simulations up to length N = 128. On the 
other hand by restricting interest to sub-manifolds of the parameter space longer chains can 
be analyzed. We have performed simulations along various lines and at points in the phase 
diagram usingwalks up to length = 1024. The schematic phase diagram conjectured by 



Kumar et al. [6 



lol 



is shown in Fig. [9] along with special lines considered in our simulations. 



To demonstrate what is estimated in a FlatPERM simulation consider for a moment a 
general polymer model with microscopic energies —£i, —£2, etc associated with configura- 
tional parameters mi, m2 respectively. Let the density of states be Ciy^rni,m2,...- Then the 
partition function is given by 

mi,m2,.-- 

where Pi = Pei, P2 = etc, and P = -j^, with ks Boltzmann's constant. FlatPERM can 
estimate CN,mi,m2,... or any sum of the CN,mi,m2,... over any number of the rrij for a range of 
lengths < Nmax- If one finds CN,mi,m2,... then one can estimate average quantities over 
this distribution for any values of Pi, P2, • • • • In our model we have mi = m and m2 = h 
with Pi = 1/T and P2 = F/T. We have performed simulations over the complete space of 
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Te 1.71 

FIG. 9: Schematic phase diagram conjectured by Kumar et al. [3, [l^]- The conjectured phase 
diagram has three different phases: 'cohapsed', 'extended' and 'stretched'. We have performed 
simulations of the whole phase space up to length Nmax = 128. The dashed lines at fixed tempera- 
tures T = 0.7 and T = 1.71 display lines along which simulations were performed for walk lengths 
up to Nmax = 1024. The six points displayed are those at which we have focused our attention. 



the variables m and h for < N^ax = 128. In this way we have estimated the density of 
states CN,m,h- We performed 10 different runs with this parametrisation for lengths up to 
Nmax = 128. We have estimated the average number of contacts per monomer {m)/N and 
the average extension per step {h)/N and their fluctuations a'^{'m)/N and a'^{h)/N. 



Previous work (see 



ll| and references therein) has estimated the ^-point to be around 



T = 1/0.663 = 1.51. With this in mind we have also performed one-parameter simulations 
with T fixed at T = 0.7 and at T = 1.71. The temperatures chosen were to ensure that one 
temperature was below and one was above the 6'-temperature as shown in Fig. [9l We also 
studied the temperature T = 2.0 with F = as a high temperature point. 

In order to delineate the possible phases we considered the points F = 0.0,0.5,1.5 for 
T = 0.7 and 1.71 and also the point T = 2.0, F = 0. In particular, we analyzed the scaling 
of the end-to-end distance R% which gives an estimate of the exponent z/. Let us start with 
F = 0. For T = 2.0 we expect that the polymer is in the extended phase with 2z/ = 1.5 and 
in Fig. [10] we find precisely that. For T = 0.7 we expect the polymer to be in the collapsed 
phase with 2z/ = 1.0 and once again our data in Fig. [TT] confirms this expectation. Now let 
us move to F = 0.5. For the low temperature T = 0.7 the series data places this point in 
the collapsed phase and the data in Fig. [11] bears this out. 



12 



However, for the point F = 0.5 with T = 1.71 the conjectured phase diagram of Kumar 
et al. {gI predicts this point to be in the extended phase with 2z/ = 1.5 while we find that 
2u = 2.00(2), so this point is in the stretched phase. For F = 1.5 and for T = 0.7 and 
T = 1.1 the conjectured phase diagram predicts a stretched phase with a value of 2z/ = 2 
and we confirm this as seen from Fig. [T2l 




0.001 0.002 0.003 0.004 0.005 0.006 
1/N 



FIG. 10: End-to-end distance divided by N^" against 1/A^ for point, {T,F) = (2.0,0.0). We see 
= 1.5(2). This point is in the extended phase. 
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FIG. 11: End-to-end distance divided by N'^'^ against N for two points, {T,F) = (0.7,0.0) and 
{T,F) = (0.7,0.6) for lengths up to = 1024. A clear linear dependence is seen implying 2i^ ~ 1. 
These points are in the collapsed phase and this exponent result is consistent with this assumption. 



It is clear that the assumption of an extended phase for F > is incorrect. However, 
while mistakenly named perhaps three phases still exist for F > 0. Considering the series 
results the other possibility is that the extended phase is indeed stretched with z/ = 1 for 
F > and that the "stretched" phase described in Kumar et al. {g], 1^ is reahy a "fully 
stretched" phase where in addition to = 1 the average height per step converges to unity: 



li,n W 



1. 



(3) 
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FIG. 12: End-to-end distance divided by N'^" against 1/N for three points, {T,F) = (0.7,1.5), 
(1.71,1.5), and (1.71,0.6). for length up to = 1024. We conclude, that all these points belong 
to a stretched phase of same type, where = 1.0. 

That is, the configurations of the polymer are essentially rod-like (with sub- dominant fiuc- 
tuations). For such rod-like configurations one would also expect in this phase that 



(m) 



lim 



0. 



(4) 



A revised conjectured phase diagram is drawn in Fig. [131 along with our lines of longer 
length simulations and the points at which we have focused our analysis. The series data 
in Fig. [7] and [8] for the low temperature regions when F = 1.2 display behavior resembling 
that delineated above for a possible "fully stretched" phase. 

To search for possible phase transitions we have estimated the fluctuations in the number 
of contacts and fluctuations in the height. For flxed force F = 1.4 the plot of the fluctuation 
against temperature in Fig. [T3] shows no sign of a growing singularity for lengths up to 

= 128 as seen in the series data for shorter lengths and smaller forces. Now we consider 
the flxed temperature lines at T = 1.71 and T = 0.7. For T = 1.71 the only sign of a 
singularity appears near F = (see Fig. [T5l) . that is the expected sign of the transition from 
the extended phase at F = to the stretched phase at F > 0. For T = 0.7 again there is 
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FIG. 13: Here is a second possible phase digram assuming there are two phases for non-zero forces, 
both with V = \. Instead of extended and stretched we have "fully stretched" and stretched. In 
such a hypothesized "fully stretched" phase the polymer is effectively in a rod-like conformation 
where the average height is approximately equal to the length of the polymer. 
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FIG. 14: The fluctuations in the number of contacts m for i*" = 1.4 fixed plotted against temper- 
ature T. There is no sign of any growing singularities which would indicate a phase transition in 
the thermodynamic limit. 

only a sign of a single phase transition in either the fluctuations of m and h (see Fig. [T5l) . 

Now the question may be asked about the nature of the peaks in the fluctuations seen 
at low temperatures at forces just above F = 1. In Fig. [T7] we plot the maximum in the 
fluctuations at fixed force for various values between F = 1.1 and F = 1.2. We note that 
while these peaks do exist they are not indicative of any divergences. Of course there may 
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FIG. 15: Plots of the fluctuations against force F for the high temperature T = 1.71. 




FIG. 16: Plots of the fluctuations against force F for the low temperature T = 0.7 

still be a weak phase transition. We now turn our attention to any possible difference in the 
conformations of the polymer in the regions labeled "stretched" and "fully stretched" . 

To test the hypothesis on which the revised conjectured phase diagram (Fig. [T^ rests 
we consider the scaling of the average height of the last monomer. In Fig. [18] we plot the 
height of the last monomer at six different points for temperatures T = 0.7 and T = 1.71. 
We observe that at the three points (T, F) = (0.7, 1.5), (1.71, 1.5) and (1.71, 0.6) the average 
height converges to a non-zero, and importantly, non-unity value. Also, at the remaining 
three points, while there are clear non-linear corrections to scaling, the average converges 
to zero. In other words no indication of a fully stretched phase can be found. A further 
test of the hypotheses leading to the revised conjectured phase diagram (Fig. [131) can be 
carried out. We assumed that for very low temperatures and large finite forces the average 
number of contacts per step goes to zero. To test this we have plotted (m) against 1/A^ for 
F = 1.2 with T = 0.1 in Fig. [191 While small (of the order of 10~®) this quantity is strictly 
increasing with length and clearly converges to a (small) non-zero value. 

nn 

We therefore conclude that the upper phase boundary proposed in [6|, [10|] does not exist 
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FIG. 17: The maximum of the fluctuations in the number of contacts m for lengths from = 32 
up to = 128 for four different forces in the region where the original series data analysis detected 
a phase transition. 
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FIG. 18: The height of the last monomer at six different points depicted in Fig. [9l We see that in 
all cases the height of the last monomers does not converge to 1, which we would expect for a fully 
stretched phase for figures at F = 1.5, T = 0.7 and T = 1.71 and for F = 0.6, T = 1.71. 



in the thermodynamic limit and the revised phase diagram in the thermodynamic limit is 
shown in Fig. [20l 
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FIG. 19: The behavior of the average number of bulk contacts {m) for length = 32 up to 
N = 128 for the point T = 0.1, F = 1.2 as a function of 1/A^. We observe an increase in the 
average number of contacts (m) with even for very low temperatures. 
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FIG. 20: The finite-size phase diagram of a self-interacting self-avoiding walk under tension in two 
dimensions for length N = 128. We have estimated the position of the maximum force Fp, and 
the force Fc for T = 0.0. The ~ 1.47. One would expect that Fc = 1.0 in the thermodynamic 
limit. 

Finally we attempt to measure the exponents associated with the collapse to stretched 
phase transition. This seems to be a second order phase transition with divergent specific 
heat. In Fig. HT] we plot the logarithm of the fluctuations in the number of contacts per 
monomer m/N against log(A^). The data in this plot is obtained at T = 0.074 and at the 
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force F for which the fluctuations are maximal. From the date we obtain estimates of the 
specific heat exponent a = 0.62(10) and the crossover exponent (p = 0.72(6). The divergence 
of the finite size specific heat is expected to be controlled by the exponent a0 and the two 
exponents are expected to by related via the scaling relation 2 — a = 1/0. 
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FIG. 21: The logarithm of the fluctuation in the number of contacts per monomer m/N against 
log(A^) at T = 0.074 for the force F at which the fluctuations are maximal. From this curve we 
can estimate the exponents a = 0.62(10) and (p = 0.72(6), which we note are not those of the 
two-dimensional 0-transition. 



V. SUMMARY 



In summary we have shown that for a model of self-interacting polymers pulled away 
from a surface in two dimension there are only two different phases for non-zero forces in 
the thermodynamic (infinite length) limit. We therefore conjecture a generic phase diagram 
as in Fig. [221 One of the phases is the collapsed phase, which is driven by the temperature 
at small forces. The other is a single stretched phase which occurs whenever the force is 
applied for temperatures higher than the ^-temperature, and for large enough forces for 
small temperatures. Importantly, the polymer is only in a fully stretched state at zero 
temperature for forces F > Fc = 1 or when the applied force is infinite. 
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FIG. 22: The final conjectured phase diagram for a two-dimensional self-interacting polymer under 
a stretching force arising from our analysis in this work. 
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